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Abstract 

New exact vacuum solutions with various singularities in the plane- 
symmetric spacetime are shown, and they are applied to the analysis of in- 
homogeneous cosmological models and colliding gravitational waves. One of 
the singularities can be true null singularities, whose existence was locally 
clarified by Ori. These solutions may be interesting from the viewpoint of the 
variety of cosmological singularities and the instability problem of Cauchy 
horizons inside black holes. 
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I. INTRODUCTION 



The important behavior of colliding gravitational waves was first shown in Khan and 
Penrose's pioneering work [1] in a plane-symmetric spacetime. The formulation for dynamics 
of colliding plane waves was given by Szekeres [2] in an extended form with non-parallel 
polarization, and the solution, which is impulsive similarly to Khan and Penrose's one, was 
derived by Nuktu and Halil [3]. 

The plane-symmetric inhomogeneous cosmological models were, on the other hand, stud- 
ied by Gowdy [4] in a vacuum spacetime with compact invariant hypersurfaces, and models 
containing space-like and time-like singularities alternately were discussed by the present 
author [5]. 

Studies on impulsive colliding gravitational waves were developed by Chandrasekhar and 
his collaborators [6] who treated not only vacuum cases but also the cases with fluids and 
electromagnetic fields, and derived also the solutions with a time-like singularity. Yurtsever 
[7] derived solutions which tend to the Kasner solutions asymptotically, contain the Cauchy 
horizons, and are connected with the Schwarzschild and Weyl vacuum solutions. 

The Cauchy horizons in the plane-symmetric spacetime were paid attention in similarity 
to the Cauchy horizons appearing in the neighborhood of the time-like singularity of the 
Ressner-Nordstrom black hole spacetime. Around the latter Cauchy horizon there is the 
extraordinary behavior of the spacetime such as the mass inflation [8] and infinite blue 
shift [9]. Recently Ori discussed the instability of the Cauchy horizon in plane-symmetric 
spacetimes in connection with that in the black hole solutions [10], and showed the possibility 
of a true null singularity in colliding gravitational waves [11]. 

In this paper we show new exact solutions representing plane-symmetric spacetimes with 
true null singularities, which are deeply related to Ori's local proof [11]. In §2, we show the 
basic equations to be solved, for the plane-symmetric spacetime s in the case of parallel 
polarization, and derive their new solutions with a kind of self-similarity. It is found by 
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analyzing one of the curvature invariants that in addition to the space-like or time-like 
singularity there appear true null singularities for a limited range of the model parameter. 
In §3, schemes of colliding waves are shown, corresponding to the above solutions, and in §4 
concluding remarks are given. In Appendices A and B the tetrad components of curvature 
tensors and some basic formulas for hypergeometric functions are shown. 



II. INHOMOGENEOUS MODEL 



The line-element of plane-symmetric spacetimes with parallel polarization is given in the 
Rosen form 



ds 2 = -2 e 2a dvdu + R [e 27 dy 2 + e" 27 dz 2 



(2.1) 



where a, R and 7 are functions of v and it, and x° = v , x 1 = u, x 2 = y and x 3 = z. The 
vacuum Einstein equations reduce to 



R,oi — 0, 



(2.2) 



27,01 + (7,0-R.i + ZiR,o)/R = 0, 



(2.3) 



where 



e 2a = iRoR^R-^e 21 , 



(2.4) 



dv( 7fi ) 2 /R fi + du( ltl ) 2 /R,_ 



R 



(2.5) 



with i?o = dR/dx°, i?i = dR/dx 1 and so on. The derivation of the above equations are 
shown in Appendix A, together with the components of the curvature tensor. 
From Eq.( |2.2| ) we have 



R = f{v)+g{u) 



(2.6) 
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where / and g are arbitrary functions. If f'iy) > and g'(u) > 0, we can put f{y) = v and 
g(u) = u without loss of generality. In the same way, if f'iy) > and g'{u) < 0, we can put 
f(v) = v and g{u) = —u. 

A. The case when the hypersurfaces R = const are space-like 

In this case we can assume f'iy) > and g'{u) > 0, and so R = v + u. Then Eqs. 
and (12.51) lead to 



27,01 + (7,Q + 7,i)/(« + u) = 0, (2.7) 

and 

[dv( lfi ) 2 + du( lA ) 2 ](v + u). (2.8) 

Here let us assume the following form of 7 : 

1 =\v\ a \u\- 1/2 A{y/u). (2.9) 

This 7 has a kind of self-similarity, because A is a function of only v/u. The parameter a 
is equal to 1/n in Ori's paper [11]. If we substitute Eq.(|2.9j) into Eq.(|2.7|), we obtain an 
ordinary differential equation for A(p) 

d 2 A dA 1 

pip + + [(2 + a)p + 1 + a}— + -(1 + 2a)A = 0. (2.10) 

This is a hypergeometric differential equation. Its independent solutions are expressed using 
hypergeometric functions as follows and have singular points at p — 0, — 1, 00 : 
For the non-integer a, we have 

A = F(l/2, 1/2 + a, 1 + a; -p) (2.11) 

and 

A = |p|- a F(l/2 - a, 1/2, 1 - a; -p). (2.12) 
For the integer a = 0, 1, 2, • • •, 
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and the second solution ( |2.12| ) is a sum of the following two solutions: 

7 oc \v \ a - 1/2 F(l/2 - a, 1/2, 1 - a; -l/p) 



(2.22) 



and 



7 oc \u\ \v\ 



-1/2 



F(l/2,l/2 + a,l + a;-l/p), 



(2.23) 



This means that the solutions are symmetric for the conversion of v — > u and u — > v. 
Here let us examine the singular behavior in the limit of v — ► 0. 

(1) Non-integer a 



For the first solution ( 2.1 1|) or ( 2.17|) in the case of a < 1/2, we find (— I) oc \v 



2a-l 



oc 



v\ 6a 3 \u\ 4 — > oo, which comes from the product (7,o) 3 ■ (7,i) 3 i n 



from Eq. (|2.15 ) and $ oc 
Eq.(Al9), so that we have the invariant curvature 7Z — > oo. If a = 1/2, we find that / and 
$ are finite, so that 1Z also is finite, and if a > 1/2, / oc \v | 2q_1 — > and $ is finite, so that 
1Z is finite. 

The second solution ( |2 . 1 2| ) or ( |2.18| ) is regular at v ~ and so / and 1Z are finite in the 
limit of v — > 0. 

It is concluded therefore that only the first solution given by fl2.11p or ( |2.17| ) for a < 1/2 
has the true null singularity. The singular behavior in the limit of u — > is the same as that 
in the limit of v — ► 0. 

As was shown in Eqs. ( |2.20|) ~ (|2.23|) , the first and second solutions at v ~ correspond 
to the first and second ones at u ~ 0, but this not the one-to-one correspondence. The 
second solution at v ~ corresponds to both solutions at u ~ 0, and the second solution at 
u — corresponds to both solutions at v ~ 0. 

As for the individual tetrad components, there are divergent components even for a > 
1/2. For example, -R(0202) and -R(0303) diverge in the limit of v — > for a < 2, because they 
include 7 j0 o oc |f| Q ~ 2 . This divergence is a kind of weak singularities. 
(2) Integer a = 0,1,2,-- • 

If a = 0, we find that 70 is finite, so that (— /), <3>, and 1Z from the first solution fl2.13p are 
finite. In the second solution (|2.14|) , the term In \p\ gives the infinity to (—1), and 7Z in the 



limit of v — > and w — > 0. If a > 1, on the other hand, 70 is finite, in both solutions (|2.13|) 
and (p. 14 ), so that TZ remains to be finite. The true null singularities appear, therefore, only 



for a = 0, if a is integer. 

Next let us consider the behavior of the solutions at R ~ 0. If we put q = p + 1, then 
Eq. ( |2.10| ) reduces to 

q(q - 1)^4 + [(2 + a)q - 1]^ + ~(1 + 2a)A = 0, (2.24) 

and the following two solutions are obtained : 

A = F(l/2, 1/2 + a, 1; g) (2.25) 

and 

A = F(l/2, 1/2 + a, 1; q) In |g| + 2^(1/2, 1/2 + a, 1; g). (2.26) 



At q ~ 0, we have F = 1 + 0(g) and Fi = 0(g), so that the first solution ( |2.25| ) is finite and 



the second ( |2.26 ) diverges as In |g| — > —00 for g| — > 0. At i? oc g ~ with finite t> and u, the 



first solution gives finite 7 and /, and so e 2a oc 1//2 from Eqs. (|2.4|) , ( 2.15|) and ( 2.16 ). If 



we define Q and T as Q = v — u and g?T = e a dR (or i? oc T 4 / 3 ), we obtain 

- 2e 2a dvdu = e 2a {-dR 2 + c/Q 2 ) oc - dT 2 + T~ 2/3 dQ 2 (2.27) 

in the neighborhood of R — 0. Accordingly the solution reduces to the Kasner solution 
with the Kasner parameter (—1/3,2/3,2/3) in the limit of R — > 0. The invariant curvature 
diverges as7^.~l/i?^oo (cf. Eq.(A19)). 

The second solution ( |2.26| ) may have a dominat contribution because of the large term 
ln|g|. It gives 7 ~ /xlni? and so e 7 oc R^, where /i = |f| a_1 / 2 . From Eq. ( p.l5j ) we obtain 
/ ~ /i 2 In |g| and so e a oc i?~ 1,/4 e / oc i?^ 2-1 / 4 . In this case, T defined by dT = e a dR gives the 
relation T oc R^ 2+3 ^, so that we have 

e a oc T Cl , i? 1/2 e 7 oc T^ 2 , i? 1/2 e" 7 oc T C3 , (2.28) 

where 



Ci = (/i 2 - l/4)/(/i 2 + 3/4), C2 = (1/2 + v)/(fi 2 + 3/4), Cs = (1/2 - /i)/(/i 2 + 3/4). (2.29) 
These exponents satisfy the Kasner relation 

E0 = E(0) 2 = 1- (2-30) 

i=l i=l 

Accordingly the solution is found to be a generalized Kasner solution at R ~ with the 
Kasner parameters depending on a spatial variable as /i = |i> | Q ~V 2 ~ |Q/2| Q_1 / 2 , where 
Q — 2t> at R ~ 0. The singularities at R = are space-like ones. 

In the limit of R — > 00 and |u| — > 00 (with v/u = const), the two solutions Q2.11|) 



and fl2.14|) satisfy the relation 7 oc R a 1//2 , and so we find for a <, =, and > 1/2 that 7 



vanishes, is finite and diverges, respectively. For a < 1/2, / in Eq. (|2.15| ) vanishes and we 



have the Kasner behavior of (—1/3, 2/3, 2/3), as in the case R — > 0. For a = 1/2, 7 and / 
remain to be constant and we have the Kasner behavior of (—1/3, 2/3, 2/3). For a > 1/2, 7 
diverges but / diverges more strongly as / oc i? 2ct_1 , so that the spacetime reduces to the 
Minkowskian spacetime with the Kasner parameter (1,0,0). 

The above treatment is applicable to the half-plane R = v + u > 0. To the half-plane 
v + u < 0, we transform v and u as v — > —v and u — > —u. Then we have R = — (v + u) > 0, 
and the same conclusion is derived in the same way. 

B. The case when the hypersurfaces R = const have time-like sections 

Now let us consider the case R = v — u. Then Eqs. Q2.3|) and (|2.5| ) lead to 



27,01 + (-7,0 + %i)/{v-u)=0, (2.31) 

and 

I = J [^( 7 ,o) 2 - ^(7,i) 2 ] {v-u). (2.32) 
If we assume 7 in the form of 

7 = \v\ a \u\- 1/2 B{v/u), (2.33) 
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we obtain an ordinary differential equation B(p): 

p(p ~ 1)^ + [(2 + a)p - 1 - a}^ + ^(1 + 2a)B = 0. (2.34) 



This equation reduces to Eq. ( |2.10| ) by putting B(p) = A(—p), so that the solutions are 



given by Eqs. ( [2.11 ) ~ ( |2.14| ) in which p is transformed to —p. It is shown similarly to 



the case A that the null singularities appear at v — and u — 0. They are true only for 
non- integer a < 1/2 and integer a = 0. 

At R ~ 0, the equation for B with respect to q = p — 1 leads to 

q(q + 1)^ + [(2 + a)g + 1]^ + ±(1 + 2a)B = 0. (2.35) 

This equation reduces to Eq. (|2.24j) by putting B(q) = A(—q), so that the solutions are 
given by Eqs. ( |2.25| ) and ( |2.26|) in which q is transformed to —q. In the limit of R — > 0, we 
find from the first solution that 7 is finite and e a oc R' 1 ^ 2 , and so we obtain 

- 2e 2a dvdu = e 2a {-dQ 2 + dR 2 ) oc - X^dQ 2 + dX 2 (2.36) 

in the neighborhood of R = 0, where Q = v + u and dX = e a dR (or R oc X 4//3 ). Accordingly 
the solution tends to the Kasner-type one with the Kasner parameter (—1/3,2/3,2/3), in 
which the time-like singularity appears in the limit of R — > or X — > 0, and the invariant 
curvature diverges as7?.~l/-R^oo (cf. Eq.(Al9)). 

From the second solution we obtain similarly the generalized Kasner-type solution ex- 
pressed by e 7 oc R^ and e a oc i? _1 / 4 e 7 oc i?^ 2 " 1 / 4 , where /1 = \v | a_1//2 . In this case, a spatial 
variable X defined by dX = e a dR gives the relation X oc R^ 2+3 ^, so that we have 

e a oc # 1/2 e 7 oc X Ca , R 1/2 e^ oc X ft , (2.37) 



where Ci)C2 ; and C3 are defined by Eq. (2.29|). Accordingly the hypersurface R = gives 



the Kasner-type time-like singularity and the Kasner parameter depends on Q - 2v , which 
is a time variable in the case B. 

In the limit of R — ► 00 and |u| — > 00 (withv/u = const), the two solutions (|2.11|) and 



( p,14j ), in which p was replaced by —p, satisfy the relation 7 oc R a 1//2 , and so we find for 



a <,=, and > 1/2 that 7 vanishes, is finite and diverges, respectively. For a < 1/2, I 
in Eq.( |2.15| ) vanishes and we have the Kasner behavior of (—1/3,2/3,2/3), as in the case 



In the above we considered the half-plane with R = v — u > 0. To the half-plane 
v — it < 0, we have to transform v and u as v — > —v and u — > —u, and use R = —v + u > 0. 
The conclusion about singularities does not depend on this transformation. The positions 
of singularities are shown in Figs. 1 and 2 for R = ±(v + u) and R = ±(v — u), respectively. 



Fig.l 



Fig.2 



III. COLLIDING GRAVITATIONAL WAVES 

Now, corresponding to two cases A and B in the previous sections, let us consider a 
spacetime containing colliding gravitational waves, which consists of the no- wave region (I), 
the outgoing- wave region (II), the incoming- wave region (III), and the colliding- waves region 
(IV). For the connection of these four regions at v = vq and u = uq, we use the Penrose 
prescription in which v and u in the metric are replaced by 

(v - v ) 6(v - v ) + v (3.1) 

and 

(U - Uq) 9(U - Uq) + U , (3.2) 
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where 8(x) is 1 and for x > and < 0, respectively. 

In the case A, the wave scheme is shown in Fig. 3, in which colliding waves are in the 
region IV within a triangle acd. The true space-like singularity is in the line ac, and for 
a < 1/2 the true null singularities are in two lines be and bf, in which one of these two lines 
can be weak singularity due to the second solution. 

In the case B, the scheme is shown in Fig. 4. Similarly colliding waves are in the region 
IV within a triangle acd. The time-like singularity ac and null singularities bf and be are 
only in the region IV. In the regions II and III weak singularities appear in the extended 
part of the above lines, because TZ vanishes there. 

In the region III, R, 7, and a are functions of only v, so that R i = 7^ = a ;1 = 0. For 
R = v ± m or R = -{v ± m ), Eqs. ([Tp, ( CT) , ( pT3|) , and (gJ§ are satisfied by an 
arbitrary function 7(f), and Eq. (All) gives 



R~ 1 / 4 



exp 



R{ lfl fdv/R fi 



(3.3) 



The tetrad components of curvature tensor vanish, except for 



R 



(0202) 



-R, 



(0303) 



-2a 



- ^7,o + 2i?(7,o) 3 - 7oo 



Here we take the functional form 



(3.4) 



1 =\v\ a \u Q \- 1/2 A(v/u ), 



(3.5) 



where A(p) and A(— p) are given by Eqs. ( 2.11| ) ~ ( |2.14j ), corresponding to the cases A and 
B, respectively. Then in the limit of R — > ( with finite values of v), we have e a oc i? -1 / 4 
and -R(o202) oc R^ 1 ^ 2 — > 00, and in the limit of v — > 0, we have for a 7^ 1/2 

,,2 



e oc exp 



a 



|2«-1 



2a 



(3.6) 



-R(0202) oc exp 



2a 2 



|2a-l 



2a- V 



ia-2 



and we have for a = 1/2 



(3.7) 



11 



e a oc M 1/4 , 



(3.8) 



-R(0202) OC V 



-2 



(3.9) 



where we assumed A(0) = 1. Accordingly, -R(o202)(= -R(0303)) diverges for a < 2, but since 
calR vanishes, these divergences give only weak singularities. 

In the connection at the boundaries of v — vo and u = uo, 7 is continuous, but 70 and 71 
are not continuous, because for example 7,0 = and 7^ at v > v and v < v®, respectively. 
Accordingly, #(0202) (= -R(0303)) has the delta function S(v — vo) and 5(u — u ) along the lines 
v — v and u = u , respectively. Moreover, since 7 oc \v | a |u| -1 / 2 and |-u| Q |t>| -1 / 2 and 70 and 
71 are not continuous at v — and u — 0, -R(o202)(= -R(0303)) nas the delta function 5{v) and 
along the lines v — and u = 0. 



In the previous sections we derived exact solutions representing inhomogeneous vacuum 
spacetimes with space-like and null singularities or time-like and null singularities, and 
applied them to the problem of colliding gravitational waves. While these null singularities 
may be closely related to the instability of Cauchy horizons, we must consider the stability of 
this null singularities themselves. For instance, there are questions whether they can coexist 



Fig.3 



Fig.4 



IV. CONCLUDING REMARKS 
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together with the electromagnetic fields and fluidal matter, and whether they can appear in 
colliding waves with non-parallel polarization. These problems will be discussed in separate 
papers. 

The physical meaning of the important parameter a may be a change rate of the shear, 
because it is included as a power index in the expression of the Kasner parameter. 

The above spacetimes are cosmologically interesting in that near the space-like or time- 
like singularity they include strong gravitational waves and at the same time show exactly the 
velocity-dominated [13] and anti-Newtonian behavior [14], in which the Kasner parameters 
depend on a spatial or time variable, respectively. 
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APPENDIX A: TETRAD COMPONENTS OF THE CURVATURE TENSOR AND 

THE RICCI TENSOR 

The following tetrads are used : 

efo = 5^, efo = tf e -, ef 2) = 5^R-^e~\ ef s) = ^iT^V (Al) 
which satisfy the condition 

9^ e ( a ) e (/3) = (a/3, (A2) 

where the non-zero components of (Ca/?) are Coi = Cio = — 1 an d C22 = C33 — 1- 

The tetrad components of curvature tensor are given by Ra^i = e ^a) e (f3) e (j) e (6)^-^^ 
and their components are expressed as 

iVoi) = -2 a i e- 2a , (A3) 
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R 



(0202) 



-2a 



(^Roo/R + 7,oo) + a (R /R + 2 7)0 ) + \(R,o/R) 2 ~ %oR,o/R ~ (7,o) 2 



(A4) 



R 



(0303) 



-2a 



(lR,oo/R - 7,oo) + a, (i2 l0 /i2 - 2 7 ,o) + -ARfi/Rf + l,oR,o/R - (l,of 



(AS) 



-R(02i2) = e 2a 



(Ir,oi/R + 7,01) + \r,oR,i/R 2 - ki,oR,i + 7,1^0) /R - 7071] , (A6) 



-R(03i3) = e 2a 



(i^oi/i? - 7,oi) + \R fi R,i/R 2 + ki,oR,i + i,iR,o)/R - 7071] , (A7) 



R 



(2323) 



-2a 



R fi R,i/R 2 - 4 7 ,o7i 



(A8) 



where _R = dR/dx°, R j01 = d 2 R/dx°dx 1 , etc. The components R(ni2) , -^(1313) and -R(m3) 
are derived from -R(o202), -R(0303) and -R(o203), respectively, by replacing d/dx° and d/dx 1 each 
other, as R >0 — > i?i and — > i? - Other components such as -R(oio2) , -R(om) etc vanish. 

The tetrad components of the Ricci tensor are given by Rt 7 s) = ( alS 'Rtayps) and their 
non-zero components are 



2(^(22) + R m ) =e- 2a R fl i/R = 0, 



(A9) 



\{R(22) - R m ) = e- 2a [2 7 ,oi + (70^,1 + 71^,0)/^] = 0, 



(A10) 



R 



(00) 



-2a 



R fi0 /R + 2a fi R, /R + -(R, /R) 2 - 2( 7 , ) 2 ] = 0, 



(All) 



R 



(ii) 



-2a 



Ru/R + 2a i-R 1/ R + ~{RJR) 2 - 2( 7 ,i) 2 ] = 0, 



(A12) 



R 



(01) 



-2a 



R 01 /R - 2a 01 + -(RM/R 2 - 27,071 



(AlS) 



From Eqs.(All) and (A12), the equations for a and ai derived and it is shown that 
their differentiations with respect to x 1 and x°, respectively give the same expression for a i 
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consistent with Eq. (A13), when we use Eqs. (A9) and (A10). Accordingly Eqs.(All) and 
(A12) are integrable for a, as in Eq. fl2.5Q in the text. 

If we use the Einstein equations (A9) ~ (A13) and the conditions R$ = const and 
R i = const, we obtain 



R 



(0101) 



"#(2323) — e 2a 



-# o#,i/# 2 + 27,07,1 



(A14) 



#(0202) — —#(0303) 



e~ 2a (R, ) 2 [ - — 7i o/ # + 2R( %0 /R fi Y - 7,oo/(# o) 2 J , (A15) 



#(1212) = -%313) = e- 2a (R^) 2 [ - ^7,i/#,i + 2#(7,i/#, - 7,n/(# 



(A16) 



1? 



(0212) — —#(0313) 



-2a 



-# o#,i/# 2 + 27071 



(A17) 



One of curvature invariants R a ^ 5 R a p^s is proportional to 



R — #(0202) #(1212) + (#(0212)) 2 + (#(0101)) 



(A18) 



Using Eqs.(Al4) ~ (A16) and Eq.(^4|) we obtain 



K = R e" 4/ $ 



(A19) 



with 



7oo/(#o) 2 - 2 J R(7o/# ) 3 + 2^7,o/#o] [7,n/(#, i) 2 - 2R( %1 /R A ) 3 



+ ^7,i/#,i 



+ 2 



1 , p/ 7,07,1 yy r 1 2 / 7,o7,i \ 2 
2i? 2+ J + L 2 i? 2+ ^#oV 



(A20) 



APPENDIX B: HYPERGEOMETRIC FUNCTIONS 



Formula ( [2.19 ) for the hypergeometric functions and their following expressions in series 



are found in the texts on special functions [12]. 

rfr, b,c;z) = i + f, a{a + 1) ' ' > + * u + J ' • ; f + * - V , (Bl) 

^ k! c(c + 1) • • • (c + K — 1) 
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F 1 (a,b,l;z) = J2 



a(a + 1) • • • (a + k - 1)6(6 + 1) • • • (b + k - 1) 



k=l 

k-1 

E 

1=0 



(k\f 



1 1 

+ 



a+l b+l l+l 



z) = (-l) c ( c -l)! z 

c-2 

*E 



(-l) fc (c- jfc-2)! 2* 



fc=0 

oo 



(a - l)(a - 2) • • • (a - c + A; + 1)(6 - 1)(6 - 2) • • • (b - c + k + 1) 



+ E 



a(a + 1) • • • (a + k - 1)6(6 + 1) • • • (6 + k - 1) 



x 



fc! c(c+ 1) • • • (c + fc- 1) 
L^v ffl + /^6 + z c + i i + iJ 



k=l 
k-1 



1=0 
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FIGURES 

FIG. 1. Inhomogeneous vacuum models with a space-like singularity at R = and null singu- 
larities at v = u = 0. The null singularities are true only for a < 1/2. Two figures are shown in 
case when (a) R = v + u and (b) R= —(v + u). 

FIG. 2. Inhomogeneous vacuum models with a time- like singularity at R = and null singu- 
larities at v = u = 0. The null singularities are true only for a < 1/2. Two figures are shown in 
case when (a) R = v — u and (b) R = —v + u. 

FIG. 3. Colliding gravitational waves with a space-like singularity and null singularities. In 
regions l(v < vq and u < no), ll(v < vq and u > no), lll(v > vq and n < no), and IV(t> > vq 
and n > no), there are the Minkowskian spacetime, the outgoing wave, the incoming wave, and 
colliding waves, respectively. These singularities are true only in the region IV. 

FIG. 4. Colliding gravitational waves with a time-like singularity and null singularities. These 
singularities are true only in the region IV in which v > v o and n > no . (a) Null singularities are in 
the lines be and bf. (b) The limiting case when one of null singularities (bf) and a null boundary 
(ad) are overlapped. 
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R = 




R = v + u 



R = - (v + u) 
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